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COT, OR COMPENSATION FOR TEMPERATURE AND ELECTRONIC GAINS 

Field of the Invention 

The present invention is directed to a method of analyzing at least two 
5 analytes at the same time using at least two fluorescent detecting entities. More 
particularly, the present invention relates to the use of fluorescently labeled 
hybridization probes to identify the genotypes at more than one nucleic acid locus by 
correcting for temperature-dependent spectral overlap of the fluorescent probes. 

10 Background and Summary of the Invention 

The continued discovery of novel genes provides a resource of genetic 
material for studying the association between genotype and disease. See Kononen, J., 
et al. (1998) Nat. Med. 4, 844-847. The majority of genetic diseases are due to single 
base alterations that may be found at multiple sites within one or several genes. See 
15 Cooper, D. N., and Krawczak, M. (1990) Hum. Genet. 85, 55-74; Neufeld, E. J. 

(1998) Hematol. Oncol. Clin. North Am. 12, 1 193-1209. For this reason, techniques 
for nucleic acid analysis often require analysis of multiple loci and sequence variants. 
For convenience, it is preferable to perform this analysis in a single reaction. 

Different colored fluorescent dyes are frequently used to increase the 
20 quantity of information obtained from a nucleic acid sample. Mansfield, E. S., et al. 
(1997) J. Chromatogr. A. 781, 295-305; Samiotaki, M., et al. (1997) Anal. Biochem. 
253, 156-161 . These dyes can be attached to primers or probes and the products 
analyzed either during PCR amplification or by post-amplification detection methods. 
Pritham, G. H., and Wittwer, C. T. (1998) J. Clin. Lig. Assay. 21, 1-9. Although post- 
25 amplification analysis increases assay time, post-amplification analysis provides a 
second level of investigation that multiplies the power of the assay. For example, 
multicolor fluorescence detection in combination with product sizing has been 
successful for identity typing using short tandem repeats, genotyping single base 
changes by minisequencing, genotyping by competitive priming. Mansfield, E. S., et 
30 al. (1997) J. Chromatogr. A. 781, 295-305; Pritham, G. H., and Wittwer, C. T. (1998) 
J. Clin. Lig. Assay. 21, 1-9; Pastinen, T., et al. (1996) Clin. Chem. 42, 1391-1397. 
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Hybridization probes provide an elegant system for homogenous PCR 
amplification and genotyping. See Lay, M. J., and Wittwer, C. T. (1997) Clin. Chem. 
43, 2262-2267; Bernard, P. S., et al. (1998) Anal. Biochem. 255, 101-107; Bernard, 
P. S., et al. (1998) Am. J. Pathol. 153, 1055-1061. In one such system, two 
5 oligonucleotide probes that hybridize to adjacent regions of a DNA sequence are 

employed, wherein each oligonucleotide probe is labeled with a respective member of 
a fluorescent energy transfer pair. In this system, a donor fluorescent dye is excited 
and transfers energy to an acceptor fluorescent dye if the two dyes are in close 
proximity, as they would be when hybridized to adjacent regions of the DNA 
10 sequence. Post-amplification melting curve analysis allows at least two alleles to be 
identified using this single pair hybridization probes. This is because variations in 
DNA sequences create mismatches with the probe resulting in characteristic shifts 
that are measured by monitoring changes in fluorescence resonance energy transfer 
during slow heating. 

15 Probes having different T m s with a variety of alleles can be multiplexed 

together to increase the power of this technique (Bernard, P. S., et al. (1998) Am. J. 
Pathol. 153, 1055-1061); however, multiplexing is ultimately limited by the number 
of T m s that can be differentiated over a range of probe melting temperatures. The 
present invention is directed to expanding the power of hybridization probe 

20 multiplexing by using color as well as T m for genotyping different alleles at one time. 
Multiple probes having different acceptor dyes are employed. Because the acceptor 
dyes emit at different wavelengths, such multicolor analysis can be used to expand the 
number of alleles which can be studied at one time. 

While multicolor analysis can be used to expand the amount of 

25 information obtained in a single melting curve, multicolor analysis requires the use of 
crosstalk compensation techniques to correct for fluorescence overlap between 
channels. These methods were originally developed for multiparameter fluorescent 
monitoring of cells using flow cytometry. Bagwell, C. B., and Adams, E. G. (1993) 
Ann. N.Y. Acad. Sci. 677, 167-184. While these techniques have proven to be useful, 

30 the temperature remains constant in flow cytometry and these methods do not correct 
for changes in temperature-dependent crossover effects. With melting curve analysis, 
the temperature ranges from 40 to 95 ° C, and significant errors can arise if 
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temperature-dependent crossover effects are not corrected. In order to adapt 
multicolor analysis for use with hybridization probes, algorithm allowances for 
different gain settings and for temperature effects on fluorescence overlap between 
channels are necessary. 
5 The present approach is to perform a calibration run during a 

temperature ramp to define the temperature dependence of each fluorophore. This 
temperature dependence is then approximated with third degree polynomials. Then, 
for fluorescent value acquisition during subsequent test runs, the temperature of the 
acquisition is used to interpolate fluorescence calibration coefficients at that 

10 temperature. Without this correction, fluorescent values would be incorrect due to 
temperature dependent spectral overlap. 

Furthermore, correction for amplifier gain is also desired, since 
experimental runs may be acquired at various gains. Because it is not practical to 
perform calibration runs for all gain combinations, the compensation file should work 

15 with experimental runs obtained at any gain setting. Correction for gain can be 

achieved by multiplying the fluorescent calibration curves by the ratio of the run gain 
to the calibration gain for each channel. 

Thus, one aspect of this invention is directed toward a method for 
determining the presence of at least two analytes using at least two fluorescent 

20 detecting entities specific for their respective analytes. The fluorescent entities are 
excited with light having the appropriate wavelength; the fluorescence is determined 
in at least two spectral channels; and the fluorescent values are corrected for spectral 
overlap and for temperature dependence of the fluorescent values. 

In another embodiment of this invention, the analytes are nucleic acid 

25 loci, and the fluorescent entities are specific for their respective loci. The fluorescent 
entities are excited with light having the appropriate wavelength; fluorescence is 
measured throughout a range of temperatures; and the signal values are compensated 
for temperature dependence. In a preferred embodiment, the fluorescent entities are 
fixed to oligonucleotides having sequences complementary to their respective loci; the 

30 oligonucleotides are annealed to the loci; and the measuring step includes monitoring 
the fluorescence during heating to melt the oligonucleotides from their respective loci. 
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The apolipoprotein ("Apo E") gene serves as a model system for the 
present invention. Single base alterations within codons 112 and 158 of the Apo E 
gene account for the three common alleles (e2, e3, and e4) and six phenotypes of Apo 
E. Mahley, R. W. (1988) Science 240,622-630. Oligonucleotide targets were 
5 synthesized to provide for adjacent hybridization probe geno typing of Apo E alleles. 
By using artificial templates, the effects of target concentration, complementary strand 
competition, probe concentration, Mg ++ concentration, and annealing conditions prior 
to melting curve analysis could be systematically studied irrespective of amplification 
technique. Thus, an additional embodiment of this invention includes optimizing 

10 target concentration and annealing conditions, and reducing complementary strand 
competition to provide optimal compensated fluorescent measurements. 

Additional features of the present invention will become apparent to 
those skilled in the art upon consideration of the following detailed description of 
preferred embodiments exemplifying the best mode of carrying out the invention as 

15 presently perceived. 

Brief Description of the Drawings 

Fig. 1 is a schematic display of fluorescent melting curve analysis for 

genotyping.; 

20 Fig. 1 A shows the equilibrium between annealing and melting at the 

T m of matched and unmatched duplexes; 

Fig. IB is a plot of time vs temperature during the melting phase of the 

reaction; 

Fig. 1C is a plot of fluorescence vs temperature for the LC Red 640 

25 probe; and 

Fig. ID is a plot of the negative derivative of the fluorescence of Fig. 

1C vs time. 

Fig. 2 shows the effect of different annealing protocols on the peak 
symmetry of derivative melting curves for a sample heterozygous at codon 158 of the 
30 Apo E gene; 

Fig. 2A illustrates a rapid cooling protocol (20° C/s) down to a 
temperature 8 0 C below the T m of a mismatched duplex; 



WO 01/12854 PCTYUS00/22083 
-5- 

Fig. 2B illustrates the rapid cooling protocol interrupted by 20 second 
holds at temperatures 8° C below the T m for the matched and then the mismatched 
duplex; 

Fig. 2C illustrates a slower temperature transition rate (1 0 C/s) down 
5 to a temperature 8 ° C below the T m of a mismatched duplex; 

Figs. 2D-F are graphs of the negative derivative (-dF/dT) melting 
curves (0.1 ° C/s) for the annealing protocols of Figs. 2A-C, respectively. 

Fig. 3 illustrates the optimization of target oligonucleotide and Mg ++ 
concentrations for genotyping codons 1 12 and 158 of the Apo E gene using the 
10 annealing protocol of Fig. 2C; 

Fig. 3 A is a -dF/dT graph measuring LC Red 705 (codon 1 12) at ImM 

MgCl 2 ; 

Fig. 3B shows LC Red 705 (codon 1 12) at 3mM MgCl 2 ; 

Fig. 3C shows LC Red 640 (codon 158) at 1 mM MgCl 2 ; and 
1 5 Fig. 3D shows LC Red 640 (codon 1 58) at 3mM MgCl 2 . Target 

concentrations were at 0.05 /uM ( — ) or 0.2 /uM (- -), fluorescein probes were at 0.1 
/uM, the acceptor probes were at 0.2 fiM, and competitor strands were not present. 

Fig. 4 illustrates competition during genotyping of Apo E targets in the 
presence of complementary strands, using the annealing protocol of Fig. 2C and 
20 plotted as the fluorescence derivative (-dF/dT) vs temperature; 

Fig. 4A shows the heterozygous melting peaks for the e3/e4 alleles 
(codon 1 12) and 

Fig. 4B shows the heterozygous melting peaks for the e2/e3 alleles 
(codon 158), in the presence ( — ) and absence (- -) of competitor strands. Each 

25 reaction included 0.1 /uM fluorescein-labeled probe, 0.2 /uM acceptor probe, 0.05 /uM 
heterozygous target with or without 0.05 /uM competitor, and either ImM MgCl 2 
(codon 1 12) or 3mM MgCl 2 (codon 158). 

Fig. 5 shows calibration and application of color compensation for 
spectral overlap with compensation of an extreme case of LC Red 640 spillover into 

30 the LC Red 705 channel; instrument gains were 1 , 30, and 50 for fluorescein, LC Red 
640, and LC Red 705 channels, respectively; the concentration of labeled probes were 
0.1 /uM fluorescein (-o-), 2 LC Red 640 (-X-), and 2 /uM LC Red 705 (— ) in 50 
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mM Tris, pH 8.5 (25° C), ImM MgCl, and 250 /iM/ml bovine serum albumin, and a 
sample containing buffer without probes (-a-) was included as a blank; 

Fig. 5A shows the original color calibration data; 

Fig. 5B shows the calibration data after color compensation; and 
5 Fig. 5C displays temperature correction to 50° C after color 

compensation; the compensation and correction algorithms were then applied to the 
LC Red 705 channel in a multiplex color assay, and fluorescence vs temperature 
curves are shown in Figs. 5D-F; the genotypes shown are homozygous 

64/homozygous e3 ( ), homozygous e3/homozygous e3 ( ), and homozygous 

10 e3/heterozygous e2/e3 (— ). 

Fig. 6 shows multiplex color genotyping of codons 112 and 158 
without competitor strand, using the annealing protocol of Fig. 2C; the 
uncompensated (Fig. 6A) and compensated (Fig. 6B) derivative melting curves 
(-dF/dT vs temperature) are shown for the LC Red 705 channel used for the analysis 
15 ofcodonll2; 

Fig. 6 A shows bleed over from the LC Red 640; 

Fig. 6B shows only the genotypes at codon 112; the genotypes shown 

are homozygous e4/homozygous e3 ( ), homozygous e3/homozygous e3 ( — ), and 

homozygous e 3 /heterozygous e2/c3 ( ). 

20 Fig. 7 shows the effect of a lower annealing temperature on peak area 

for the mismatched g3 in genotyping of the heterozygous e2/e3 sample; 

Fig. 7A shows genotyping after cooling at 1 ° C/s to 42° C, the 
temperature used for multiplexing in Fig. 6; and 

Fig. 7B shows genotyping after cooling at 1 ° C/s to 48 ° C, the 
25 temperature used for genotyping optimization at codon 158. 

Detailed Description of the Invention 

In describing and claiming the invention, the following terminology 
will be used in accordance with the definitions set forth below. 
30 As used herein, "nucleic acid," "DNA," and similar terms also include 

nucleic acid analogs, i.e. analogs having other than a phosphodiester backbone. For 
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example, the so-called "peptide nucleic acids," which are known in the art and have 
peptide bonds instead of phosphodiester bonds in the backbone, are considered within 
the scope of the present invention. 

As used herein, "fluorescence resonance energy transfer pair" or 
5 "FRET pair" refers to a pair of fluorophores comprising a donor fluorophore and 
acceptor fluorophore, wherein the donor fluorophore is capable of transferring 
resonance energy to the acceptor fluorophore. In other words the emission spectrum 
of the donor fluorophore overlaps the absorption spectrum of the acceptor 
fluorophore. In preferred fluorescence resonance energy transfer pairs, the absorption 
1 0 spectrum of the donor fluorophore does not substantially overlap the absorption 
spectrum of the acceptor fluorophore. 

As used herein, "a donor oligonucleotide probe" refers to an 
oligonucleotide that is labeled with a donor fluorophore of a fluorescent resonance 
energy transfer pair. 

1 5 As used herein, "an acceptor oligonucleotide probe" refers to an 

oligonucleotide that is labeled with an acceptor fluorophore of a fluorescent resonance 
energy transfer pair. 

As used herein, "FRET oligonucleotide pair" refers to the donor 
oligonucleotide probe and the acceptor oligonucleotide probe pair that form a 

20 fluorescence resonance energy transfer relationship when the donor oligonucleotide 
probe and the acceptor oligonucleotide probe are both hybridized to their 
complementary target nucleic acid sequences. 

The present invention is directed a method for using at least two 
fluorescent detecting entities for determining the presence of at least two analytes. 

25 While this technique may be applied to a variety of uses, a preferred embodiment of 
this invention is a method of screening multiple loci of nucleic acid sequences for the 
presence of mutations or polymorphisms. More particularly, the present invention 
allows for a rapid procedure that can be entirely conducted within a single reaction 
vessel, for detecting mutations and polymorphisms at multiple loci of a genomic DNA 

30 sample prepared from an individual organism. While the method may be used on any 
suitable nucleic acid sample, the method can be performed on PCR amplification 
products and can be conducted within the same single reaction vessel as used for PCR 
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amplification. The preferred method of this invention comprises the steps of 
combining a biological sample comprising nucleic acid sequences with two or more 
FRET oligonucleotide pairs, illuminating the biological sample with the appropriate 
wavelength, monitoring the fluorescence as a function of temperature during heating, 
5 and compensating the fluorescent values. 

Fluorescent probes for use in detecting and monitoring DNA include 
double-stranded-DNA-specific dyes and sequence-specific probes. Fig. 1 diagrams a 
hybridization scheme based on resonance energy transfer between fluorophores on 
two adjacent probes. As further defined in Table 1, this method is sequence specific 

10 and allows analysis with melting curves. When the two labeled probes are hybridized 
to the same template strand, the two members of the FRET pair are brought into close 
proximity and energy transfer can occur (Fig. 1A). Because of a mismatch between 
one of the alleles and one of the oligonucleotide probes, melting will occur at a lower 
temperature for that allele, thus separating the members of the FRET pair. The allele 

15 containing the mismatch will have a characteristically lower T m than the allele that 
matches the probe (Figs. 1B-D). 

Hybridization probes specific for different alleles and labeled with 
different acceptor colors may be used simultaneously to multiplex fluorescent melting 
curves for genotyping. Two separate FRET oligonucleotide pairs, each specific for 

20 one locus and each comprising a different acceptor dye may be used at the same time. 
Acceptable fluorophore pairs for use as fluorescent resonance energy transfer pairs are 
well know to those skilled in the art and include, but are not limited to, 
fluorescein/rhodamine, phycoerythrin/Cy7, fluorescein/Cy5, fluorescein/Cy5.5, 
fluorescein/LC Red 640, and fluorescein/LC Red 705. 

25 The nucleotide sample to be analyzed may be PCR amplification 

products provided using the rapid cycling techniques described in US Patent No. 
5,455,175. Rapid temperature cycling is contrasted to conventional temperature 
cycling, wherein it is shown that 30 cycles of amplification can be completed in 15 
minutes and the resulting PCR products contain many fewer side products. Thus, 

30 with rapid cycling the required times for amplification are reduced approximately 
1 0-fold, and specificity is improved. 
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Real time monitoring of melting curves during slow heating can be 
conducted through the use of fluorescent dyes and fluorescently labeled probes. The 
use of real time fluorescent monitoring of nucleic acid reactions have been described 
in U.S. Patent No. 5,455,175, issued on October 3, 1995, and U.S. Patent Applications 
5 Serial Nos. 08/869,275 and 08/869,276, each filed on June 4, 1 997. 

Fluorescent melting curve analysis is a rapid and effective method for 
genotyping. See, for example, Lay, M. J., and Wittwer, C. T. (1997) Clin. Chem. 43, 
2262-2267; Bernard, P. S., et al. (1998) Anal. Biochem. 255. 101-107; Bernard, P. S., 
et al. (1998) Am. J. Pathol. 153, 1055-1061; Lyon, E.. et al. (1998) Mol. Diag. 3, 203- 

10 210. Rapid cycle PCR and single base genotyping can be done in a single tube and 
without intervening steps within 40 minutes. The algorithms and methods of the 
present invention demonstrate that solution color multiplexing is possible with 
fluorescence-capable temperature cyclers such as the LightCycler™. See Wetmur J. 
G. (1995) in Molecular Biology and Biotechnology (Meyers, R. A., ed.), pp. 605608, 

15 VCH Publishers, Inc., New York, NY. This multiplexing technology could also be 
applied to quantification methods, such as those using internal controls for 
competitive PCR. Sidhu, M. K., et al. (1998) in Gene Quantification (Ferre, F., Ed.), 
pp. 265-276, Birkhauser, Boston, MA. 

The Apo E gene serves as a model system for the present invention. 

20 Single base alterations within codons 1 12 and 1 58 of the Apo E gene account for the 
three common alleles (e2, e3, and e4) and six phenotypes of Apo E. Mahley, R. W. 
(1988) Science 240, 622-630. Oligonucleotide targets were synthesized to provide for 
adjacent hybridization probe genotyping of Apo E alleles. The FRET pairs used are 
fluorescein/LC Red 640, and fluorescein/LC Red 705. The sequences are shown 

25 below in Table I. In addition to compensation for the temperature dependence of the 
fluorescent dyes, effects of target concentration, complementary strand competition, 
probe concentration, Mg + " concentration, and annealing conditions prior to melting 
curve analysis are studied and optimized. 
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Table I. 

Oligonucleotides and Probes Used for Genotyping the Model Apolipoprotein E Locus 
Sequences for genotyping codon 1 12 

5 e3 Target GGCGCAGGCCCGGCTGGGCGCGGACATGGAGGA CGTGTGCG 
GCCGCCTGG TGCAGT a [SEQ ID NO: 1 ] 
eA Target GGCGCAGGCCCGGCTGGGCGCGGACATGGAGGA CGTGCGC 

GGCCGCCTGG TGCAGT a [SEQ ID NO:2] 
Fluorescent CCAGGCGGCCGCACACG-fluorescein [SEQ ID NO:3] 
10 probes LCT705-CCTCCATGTCCGCGCCCAGCCGGGCCTGCG [SEQ ID 

NO:4] 



Sequences for genotyping codon 158 
e2 Target GCGGCTCCTGCCCGATGCCGATGACCTGCA GAAGTGCCTGG 

CAGTGT ACCA a [SEQ ID NO:5] 
e3 Target GCGGCTCCTGCCCGATGCCGATGACCTGCA GAAGCGCCTGG 

CAGTGT ACCA a [SEQ ID NO:6] 
Fluorescent ACACTGCCAGGCACTTC-fluorescein [SEQ ID NO: 7] 
probes LCR640-GCCAGGTCATCGGCATCGGGGCAGGAGCC [SEQ ID 

NO:8] 



"Underlined area indicates fluorescein probe target region; bold face indicates 
location of mismatch. 

25 The rate of oligonucleotide probe/target hybridization can be 

determined by a single nucleation step. However, probe/target hybridization is 
influenced by salt concentration, temperature, and secondary structure. Wetmur, J. G. 
(1995) in Molecular Biology and Biotechnology (Meyers, R. A., Ed.), pp. 605-608, 
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VCH Publishers, Inc., New York, NY. It is believed that these same conditions 
influence the quality of genotyping during probe/target melting. For instance, 
equivalent fluorescent melting peak areas of a heterozygous sample would require 
having nearly the same concentrations of probe annealed to each target prior to 
5 melting. Thus, for genotyping, annealing temperatures are selected for optimizing 
rates of probe/target hybridization. Empirically, maximum probe/target hybridization 
rates occurred 5-10° C below the duplex T m . Presently, a temperature approximately 
8° C below duplex T m is used as the target temperature. 

Using target annealing temperatures specific for each site, a gradual 

10 transition rate before melting provided the best genotyping (Fig. 2C). In addition, the 
MgCl 2 and target concentrations affect the genotyping quality at each site. A high 
MgCl : concentration produces poor peak differentiation for heterozygous genotyping 
at codon 1 12, but optimal genotyping at codon 158 in the presence of a low target 
concentration. Genotyping at codon 158 depends on a low target concentration for 

1 5 better representation of the mismatched duplex. One possible explanation for this is 
that the higher ratio of probe to target may better insure saturation of both strands of 
targets prior to melting. 

The finding that the sites optimize to different template concentrations 
suggests that the cycle at which amplification is stopped and investigation of the 

20 sequence begins is important. Melting curve analysis usually begins after the plateau 
phase of amplification is reached. See Bernard, P. S.. et al. (1998) Anal. Biochem. 
255, 101-107; Bernard, P. S., et al. (1998) Am. J. Pathol. 153, 1055-1061. However, 
this point may be too late for optimal genotyping at sites such as codon 158. As an 
alternative, melting curve analysis can be programmed to begin at any cycle to provide 

25 optimal genotyping. 

A single salt concentration and target concentration was chosen for 
multiplexing hybridization probes by color. A low salt concentration was chosen due 
to codon 112 and a low template concentration was chosen due to codon 158. The 
lower MgCl 2 concentration decreased the area of the mismatched peak at codon 158. 

30 However, this effect was partly corrected by dropping the annealing temperature for 
multiplexing farther than 8° C below the T m , the temperature initially selected. The 
lower annealing temperature used for multiplexing was not necessary for single-site 
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analysis at codon 158 since equivalent heterozygous peaks could be produced with a 
higher MgCl 2 concentration. In general, it may be necessary to use annealing 
temperatures lower than 8° C below the T m to increase probe/target annealing under 
low MgCl 2 concentrations. 
5 The plateau phase of amplification is to a large extent due to the 

accumulating complementary product strands competing for primer annealing sites. 
Wittwer, C. T, et al. (1997) BioTechniques 22, 130-138. This would be consistent 
with the observed competition during genotyping between the complementary strand 
and hybridization probes for the target. Although this competition does not prevent 

10 melting curve analysis of germ line mutations, it may limit the detection sensitivity of 
somatic mutations within a background of wild-type product. Various techniques 
have been developed to circumvent the competition from the nontarget strand during 
analysis. For example, peptide nucleic acid (PNA) probes effectively compete against 
the complementary strand by forming an extraordinarily stable duplex with the target. 

15 Egholm, M., et al. (1993) Nature 365, 566-568. Another alternative is to bind the 
target to solid support and ihen remove the competitor strand. Pastinen, T., et al. 
(1996) Clin. Chem. 42, 1391-1397. Asymmetric PCR amplification prior to melting 
curve genotyping would reduce competition while maintaining a homogenous assay. 
In addition, other techniques directly produce single stranded product, such as nucleic 

20 acid sequence-based amplification (NASBA) or rolling-circle amplification. Reitsma, 
P. H., et al. (1996) Blood Coagul. Fibrinolysis 7, 659-663; Lizardi, P. M., et al.(1998) 
Nat. Technol. 19, 225-232. 

Prior art color compensation algorithms, developed for flow cytometry, 
do not compensate for changes in signal gains or temperature-dependent crossover 

25 effects. The change in crossover coefficients with temperature is not an issue when 
the temperature remains constant. However, with melting curve analysis, the 
temperature ranges from 40-95° C and significant errors arise if the calibration data 
used to derive the crossover coefficients are obtained at a different temperature than 
the data to be evaluated. The approach of the present invention is to acquire a 

30 calibration run over a slow temperature ramp and to calculate temperature-specific 
crossover coefficients at each temperature by interpolation of a 3 rd degree polynomial. 
That is, a temperature-specific crossover matrix is calculated at the temperature of 
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each data point to be evaluated. Similarly, gain correction is possible by calculating a 
new matrix from fluorescence values that are multiplied by the gain ratio of the data 
point to the calibration run. As long as the electronic signal gains are accurate, the 
new crossover coefficient matrix is matched to the experimental data. 

Color compensation for spectral overlap of fluorescent dyes has been 
modified from flow cytometry techniques. See Bagwell, C. B., and Adams, E. G. 
(1993) Ann. N.Y. Acad. Sci. 677, 167-184. Custom analysis software was written in 
LabView (National Instruments, Austin, TX). In summary, a calibration run of each 
pure fluorescent oligonucleotide and autofluorescence control is first obtained with 
fluorescence values acquired in each channel. Signal crossover constants are 
calculated from these values and used to convert observed fluorescence (o) to actual 
signal fluorescence (s) by matrix algebra. 

Several modifications to the color compensation algorithms developed 
for flow cytometry were necessary for use on genotyping data. Following the notation 
given in Bagwell, C. B., and Adams, E. G. (1993) Ann. N.Y. Acad. Sci. 677, 167- 
1 84, hereby incorporated by reference, and eliminating the acquisition of multiple 
events (relevant to flow cytometry but not to solution fluorescence), the crossover 
constants are: 

E bo". ") - 

n = \ 

where: k(i, j) = crossover signal of dye i in channel j 

o(i, j) = observed signal of dye / in channel j 

a(j) = autofluorescence in channel j 

n = channel index 

N = maximum channel 

Preliminary experiments have shown that the fluorescence of some of 
the dyes, and the crossover constants of some dye/channel combinations, are 
temperature dependent over the desired range of analysis (40-95° C). Therefore, 
calibration runs were obtained by acquiring fluorescence values continuously during a 
0.2° C/s temperature ramp from 40 to 95° C. The temperature vs fluorescence curves 
were nearly linear, but were better approximated by 3 rd degree polynomials. The 3 rd 
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degree polynomial coefficients for the temperature vs fluorescence curves of each of 
the dye/channel combinations and the autofluorescence controls of each channel were 
stored for temperature interpolation. 

To color compensate fluorescence data, the temperature of each 
5 acquisition was used to interpolate temperature-specific fluorescence from the 

calibration curves. These values were then adjusted for any change in electronic gain 
between the calibration and data run as follows: 

o(i,j) = [w(i,j)T 3 +x(i,j)T 2 +y(i,j)T + zfi,j)] [G D (j)/G c (j)] 
a(j) = [m(j)T + nO)T 2 +p(j)T+ q(j)} [G D (j) / G c (j)} 
1 0 where: w(i, j), x(i, j), y(i, j), and z(i, j) are 3 rd degree polynomial coefficients for the 

temperature vs fluorescence curve of dye / in channel j, and m(j), n(j),p(j), and 
q(j) are 3 rd degree polynomial coefficients for the temperature vs fluorescence 
autofluorescence curve in channel j 
T = acquisition temperature 
15 G D (j) = gain of channel j during the data run 

G c (j) = gain of channel j during the calibration run 

The actual signal fluorescence of each dye is calculated by the matrix 

equation: 

S = K'[0-A] 

20 where: S = actual signal fluorescence of each dye 

K = crossover constant of each dye in each channel 
O = observed fluorescence in each channel 
A = observed autofluorescence in each channel 
If temperature correction is desired as well as color compensation, the 
25 calibration data are color compensated by the above procedure and the temperature vs 
actual signal fluorescence curves are fit to 3 rd degree polynomials and the coefficients 
stored for temperature interpolation. The data run is color compensated and the 
fluorescence of each dye during each acquisition is temperature corrected as follows: 
% = s [sJTs) s : (T :> )\ 
30 where: s TC = temperature corrected signal fluorescence 
5 = signal fluorescence 

s c (Ts) = interpolated signal fluorescence from the calibration run at a 
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standard temperature (selected by the user) 
s c (T q) = interpolated signal fluorescence from the calibration run at the 
data acquisition temperature 
An example of the results of color compensation and temperature 
5 correction is shown in Fig. 5. 

Example 1 
Preparation of the Oligonucleotides 

The human apolipoprotein E gene sequence (GenBank accession 
K00396) was used to design unlabeled target strands, unlabeled competitor 

10 (complementary) strands, and fluorescently labeled probes. The sequences are shown 
above in Table 1 . All oligonucleotides were synthesized by standard phosphoramidite 
chemistry (Pharmacia Biotech Gene Assembler Plus, Piscataway, NJ). The 5'-labeled 
LC Red 640 (LightCycler Red 640) and LC Red 705 (LightCycler Red 705) acceptor 
probes were synthesized using a C6dT amino-modifier (Glen Research, Sterling, VA) 

15 and an LC Red 705 phosphoramidite (Roche Molecular Biochemicals), respectively. 
The LC Red 640 N-hydroxysuccinimide ester (Roche Molecular Biochemicals) was 
manually conjugated to an equimolar concentration of oligonucleotide in 0.1 M borate 
buffer, pH 9.0. Incubation was performed for two hours in the dark. Both acceptor 
probes were synthesized with a chemical phosphorylation reagent (Glen Research) at 

20 the 3' end. The 3'-fluorescein-labeled probes were synthesized on fluorescein- 
controlled pore glass cassettes (BioGenex, San Ramon, CA). 

A S'-trityl group was retained during synthesis on the oligonucleotide 
target and competitor strands, the pre -conjugated LC Red 640 probe, and the 
fluorescein-labeled probes. Full-length sequences for all oligonucleotides and probes 

25 were purified by reverse-phase C 1 8 high-pressure liquid chromatography (HPLC) 
using a 3.5 p Waters Symmetry Column 4.6 x 150-mm (Waters. Milford, MA). The 
mobile phase consisted of 0.1 mol/L triethylammonium acetate, pH 7.0, and a 10 to 
30% (fluorescein probe) or 10 to 80% (unlabeled oligonucleotides and LC Red 705 
probe) gradient of acetonitrile (1 ml/minute). The eluate was monitored with tandem 

30 absorbance and fluorescence detectors (Waters 486 and 474, Milford, MA). For 
probes, fractions with coincident A 260 and fluorescence peaks were collected. Final 
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detritylation was performed on a Polypack column (Glen Research, Sterling, VA) and 
products were eluted with 50% acetonitrile. 

The oligonucleotides and probes were vacuum dried, resuspended in 1 
mM Tris-HCl, 0.1 mM EDTA. pH 8.0 and then quantified. Unlabeled 
oligonucleotides with A 260 /A 280 values between 1.6 and 1.8 were used. The 
concentration ratio of fluorophore to oligonucleotide was calculated for each probe. 
Probes with ratios outside 0.9 to 1.1 were further purified on a denaturing 15% 
polyacrylamide gel containing 7 M urea followed by HPLC to remove the urea. 

A 3 '-fluorescein labeled 17-mer probe was used at each codon for 
genotyping. Although these probes were different in sequence, both formed an A:C 
mismatch flanked by C:G pairs when hybridized to the non-complementary target 
(Table 1). The mismatch in both duplexes was positioned 5 base pairs in from the 3'- 
end of the probe. 

Genotyping at each site was optimized for reagent concentrations and 
annealing conditions prior to melting curve analysis (see Examples 2 and 3 below). 
Concentrations of MgCl 2 , target strand, competitor strand, and fluorescein probe were 
varied. The MgCl 2 concentrations were varied from 1 to 5 mM. The concentrations 
of target oligonucleotides were varied from 0.05 to 0.2 uM and competitor 
oligonucleotides were at an equimolar concentration to the target strand. Fluorescein- 
labeled probes were varied from 0.1-0.2 uM. Each 10 uL reaction also contained 50 
mM Tris, pH 8.5 (25 °C), 250 ug/ml bovine serum albumin and 0.2 uM of the 
acceptor probe (LC Red 640 or LC Red 705). Fig. 2 shows the different temperature 
conditions used for probe annealing prior to genotyping. Heterozygous samples were 
used to determine optimal conditions for genotyping, assessed as similar shaped 
Gaussian curves and equivalent peak areas. 

A single temperature protocol and common reagents were necessary 
for multiplex genotyping. The reagents used were 1 mM MgCl 2 , 0.05 uM target 
oligonucleotides for each site, 0.2 uM acceptor probe at each site, 0.1 uM fluorescein 
labeled probe spanning codon 112 and 0.2 uM fluorescein labeled probe spanning 
codon 158. The temperature protocol for multiplex genotyping is given in Fig. 2C. 

The samples were loaded into composite plastic/glass capillary 
cuvettes, capped, briefly centrifuged and loaded into a 32-sample carousel of a 
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thermal cycler with 3-color fluorescence monitoring capability (LightCycler™, Roche 
Molecular Biochemicals). Twelve samples were included in each melting curve 
analysis run. During the slow heating (0.1° C/s) phase of the genotyping protocol, 
fluorescence was continuously monitored with 100 ms acquisitions. The fluorescence 
5 values of fluorescein, LC Red 640, and LC Red 705 was acquired in exact temporal 
coincidence using a linear arrangement of dichroic bandpass filters. See Wittwer, 
C. T.,et al. (1997) BioTechniques 22, 176-181. The data were displayed as 
fluorescence (F) vs temperature (T) and as derivative fluorescence (-dF/dT) vs 
temperature plots, as shown in Figs. 1C and ID. 

10 Example 2 

Optimization of Annealing Conditions 
The quality of melting curve analysis depends not only on the melting 
protocol, but also on the annealing conditions before the melting curve acquisition. 
Fig. 2 shows that slow and stepwise temperature annealing protocols result in more 

15 symmetrical heterozygous melting peaks then annealing using a single rapid 

temperature transition. This result was observed even though the total annealing 
times were approximately the same (compare Figs. 2A. B, and C), suggesting that 
linear cooling through the T m s of each duplex favors equal formation of each duplex. 
Table 2 shows a comparison of the T m shifts created by the same 

20 mismatch at the two sites using different MgCl 2 concentrations. The mismatches 
destabilize the probe/target duplexes at each site by approximately the same amount. 
Furthermore, the AT m between the matched and mismatched duplexes was increased 
with a higher MgCl 2 concentration, while increasing the fluorescein probe 
concentration increased the T ra s of both duplexes without affecting the AT m (data not 

25 shown). 
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20 
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A:C 


64.2 
57.2 
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The clear resolution of a heterozygote at codon 1 12 (spanned by an 82% 
GC probe) was dependent on a low MgCl 2 concentration (Fig. 3 A vs 3B). In contrast, 
25 heterozygote genotyping at codon 158 (spanned by a 59% GC probe) was more symmetric 
with a higher MgCl, concentration and a low target concentration (Fig. 3C vs 3D). 

Genotyping from a double stranded product (i.e., equal concentrations of 
target and complementary strand) is illustrated in Fig. 4. In the presence of competitor 
(complementary strand), genotyping optimized to the same parameters shown in Fig. 3 
30 but with a reduction in fluorescence signal to approximately 60 % at each site. 

Example 3 

Color Compensation and Temperature Correction 

Fig. 5 displays the effect of color compensation and temperature correction 
applied to a calibration run and to color multiplexing of the Apo E locus. The 

35 temperature dependence of LC Red 705 is substantial, whereas LC Red 640 fluorescence 
is almost constant with temperature (Fig. 5A). As expected, only the LC Red 705 
calibration sample shows fluorescence after color compensation (Fig. 5B) and the 
temperature dependence is removed after temperature correction (Fig. 5C). The original 
fluorescence vs temperature traces for the LC Red 705 channel (codon 1 12) are complex 

40 and not interpretable because of fluorescence overlap from LC Red 640 (Fig. 5D). 

However, single temperature transitions are discerned after color compensation (Fig. 5E). 
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After temperature correction (Fig. 5F), the baseline slope outside of the melting transition 
region is increased. 

Multiplexing by color and T m provides simultaneous identification of 
variants at codons 1 12 and 158 in a single tube within 10 minutes. Fig. 6 shows the 
5 uncompensated (A) and compensated (B) fluorescent melting peaks for genotyping codon 
112. Without compensation, there is substantial bleed over from the LC Red 640 acceptor 
dye into the channel for monitoring the LC Red 705 acceptor dye (Fig. 6A). This results 
in apparent products with T m s of 55° C and 62° C. However these additional peaks 
disappear after compensation, revealing the true homozygous genotypes at codon 1 12. 

10 All the genotypes for the 6 naturally occurring Apo E protein isoforms could be 

distinguished using this method. In addition, the other 3 possible genotypes not normally 
found in human populations were constructed using synthesized templates and were 
correctly analyzed (data not shown). 

Multiplex genotyping was optimal using low concentrations of both target 

15 and MgCl 2 . Using the lower MgCl 2 concentration for genotyping codon 158 decreased 
the peak area of the mismatched duplex. Nevertheless, this resulting asymmetry in the 
heterozygous peaks could be partially compensated for by lowering the annealing target 
temperature (from 48° C to 42° C) before melting (Fig. 7). 

Although the invention has been described in detail with reference to 

20 preferred embodiments, variations and modifications exist within the scope and spirit of 
the invention as described and defined in the following claims. 
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CLAIMS : 



1 . A method for determining the presence of at least two analytes in a 
sample, comprising: 

5 providing at least two detecting entities, each of said entities comprising a 

fluorescent label different from the fluorescent label of any of the other detecting entities, 
and each of said entities capable of specific binding with their respective analytes, 

contacting the sample with the detecting entities to allow specific binding 
of the analytes and detecting entities, 
1 0 exciting each of said labels with light having an appropriate wavelength 

for inducing fluorescence, 

determining values of fluorescence of the labels in at least two different 
spectral channels, and 

compensating said values for spectral overlap, wherein the compensating 
15 step includes correcting for temperature dependence of the fluorescent values. 

2. The method of claim 1 wherein the determining step is performed 
throughout a range of temperatures. 

3. The method of claim 2 wherein the correction for temperature 
dependence is performed through calculating a temperature specific overlap coefficient. 

20 4. The method of claim 1 wherein the analytes are nucleic acid 

amplification products. 

5. The method of claim 4 further comprising the step of determining 
the genotype of each of the nucleic acid amplification products. 

6. The method of claim 1 further comprising the step of compensating 
25 said values for amplifier gain. 

7. A method of analyzing the genotype of multiple loci of nucleic acid 
sequences comprising the steps of: 

providing a sample of nucleic acid having a first locus and a second locus, 
providing a first detecting entity, the first detecting entity comprising a first 
30 fluorescent label and having specificity for the first locus, 

providing a second detecting entity, the second detecting entity comprising 
a second fluorescent label and having specificity for the second locus, 
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contacting the sample of nucleic acid with the first and second detecting 

entities, 

exciting the first and second entities with light having the appropriate 
wavelength to induce fluorescence of the first and second fluorescent labels, 

measuring throughout a range of temperatures a first fluorescent signal of 
the first detecting entity in a first spectral channel and a second fluorescent signal of the 
second detecting entity in a second spectral channel, and 

compensating the first and second signal values for temperature 
dependence of the fluorescent signals. 

8. The method of claim 7 wherein the compensating step is performed 
through calculating a temperature dependent coefficient. 

9. The method of claim 8 further comprising the step of compensating 
for amplifier gain. 

10. The method of claim 7 wherein the first detecting entity further 
comprises a sequence of nucleic acid complementary to at least a portion of the first locus 
and the second detecting entity further comprises a sequence of nucleic acid 
complementary to at least a portion of the second locus. 

1 1 . The method of claim 1 0 wherein the contacting step includes 
annealing the first and second detecting entities to their respective loci, and the measuring 
step includes monitoring fluorescent signals in the spectral channels during heating to 
melt the detecting entities from their respective loci. 

12. The method of claim 1 1 wherein the contacting step is performed 
by lowering the temperature by 1 ° C/s. 

1 3 . The method of claim 7 wherein the first fluorescent label comprises 
a first donor dye, a first acceptor dye, and a first pair of oligonucleotides, each 
oligonucleotide attached to one of the dyes, and each oligonucleotide substantially 
complementary to neighboring regions of the first locus. 

14. The method of claim 13 wherein the second fluorescent label 
comprises a second donor dye, a second acceptor dye, and a second pair of 
oligonucleotides, each oligonucleotide attached to one of the dyes, and each 
oligonucleotide substantially complementary to neighboring regions of the second locus. 
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15. The method of claim 14 wherein the first donor dye is fluorescein 
and the first acceptor dye is LC Red 705. 

16. The method of claim 15 wherein the second donor dye is 
fluorescein and the second acceptor dye is LC Red 640. 

17. The method of claim 13 wherein the first locus has multiple alleles 
and the Mg ++ concentration is optimized to produce relatively even melting peaks for the 
alleles when the compensated signal values are plotted as -dF/dT vs temperature. 

18. The method of claim 1 3 wherein the first locus has multiple alleles 
and the concentration of the nucleic acid comprising the first locus is optimized to 
produce relatively even melting peaks for the alleles when the compensated signal values 
are plotted as -dF/dT vs temperature. 
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